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Abstract 

Background: Human brucellosis incidence in China was divided into 3 stages, high incidence (1950-1 960s), decline 
(1970-1 980s) and re-emergence (1 990-2000s). Human brucellosis has been reported in all the 32 provinces, of which 
Inner Mongolia has the highest prevalence, accounting for over 40% of the cases in China. To investigate the 
etiology alteration of human brucellosis in Inner Mongolia, the species, biovars and genotypes of 60 Brucella isolates 
from this province were analyzed. 

Methods: Species and biovars of the Brucella strains isolated from outbreaks were determined based on classical 
identification procedures. Strains were genotyped by multi locus sequence typing (MLST). Sequences of 
9 housekeeping genes were obtained and sequence types were defined. The distribution of species, biovars 
and sequence types (STs) among the three incidence stages were analyzed and compared. 

Results: The three stages of high incidence, decline and re-emergence were predominated by B. melitensis biovar 
2 and 3, B. abortus biovar 3, and B. melitensis biovar 1, respectively, implying changes in the predominant biovars. 
Genotyping by MLST revealed a total of 14 STs. Nine STs (from ST28 to ST36), accounting for 64.3% of all the STs, were 
newly defined and different from those observed in other countries. Different STs were distributed among the three 
stages. ST8 was the most common ST in 1950-1 960s and 1990-2000s, while ST2 was the most common in 1970-1 980s. 

Conclusions: The prevalence of biovars and sequence types of Brucella strains from Inner Mongolia has changed over 
time in the three stages. Compared with those from other countries, new sequence types of Brucella strains exist 
in China. 
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Background 

Brucellosis is a zoonotic disease worldwide. It is classified 
as a reemerging infectious disease because of the increas- 
ing incidences in recent years. The genus Brucella, the 
etiological pathogen of brucellosis, is divided into at least 
ten species, including the six classical species: B. abortus 
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(bovine), B. melitensis (caprine and ovine), B. ovis (ovine), 
B. canis (canine), B. suis (porcine), and B. neotomae on the 
basis of host specificity, antigenic differences and bio- 
chemical characteristics [1,2]. The classical species are di- 
vided into biovars even though the differences defining 
some biovars are minor and somewhat subjective. Four 
out of the six classical species are pathogenic for humans: 
B. abortus, B. melitensis, B. suis and rarely B. canis. B. 
melitensis and B. abortus are highly pathogenic and a fre- 
quent cause of human brucellosis. 

Brucellosis is epidemic in China, and cases have been re- 
ported in all 32 provinces [3,4]. Since the 1950s, systematic 
surveys on the prevalent distributions and epidemiological 
characteristics of brucellosis have been carried out in China 
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[4]. Human brucellosis was highly endemic from the mid- 
1950s well into the 1970s, and then incidence decreased 
until the mid-1990s, after which it increased sharply. Based 
on this change in incidences, brucellosis in China was di- 
vided into three periods, high incidence (1950-1960s), de- 
cline (1970-1980s) and re-emergence (1990-2000s) [4]. 
Although high incidence has been reported in many prov- 
inces, the etiological pathogen has not been genetically 
characterized. The main species, biovars, and genotypes of 
the Brucella isolates and their changes in the three differ- 
ent periods remain largely unknown. Of the epidemic foci, 
the Inner Mongolia Autonomous Region had the highest 
incidence, accounting for about 40% of the reported cases 
[5]. The incidence and epidemiology of brucellosis in this 
province could represent characteristics of brucellosis in 
China [4]. Spatial and temporal distributions of human 
brucellosis in this province during the time period of 
1999-2008 were systematically investigated in a recent 
study. The annual incidence varied greatly from 0 to 
818.52/100,000 at county levels during the study period. 
The highest incidence of the disease occurred in Abaga 
County. The spatial distribution of the disease clustered in 
the northeastern and central districts [5]. 

Multi locus sequence typing (MLST) is a DNA 
sequence-based typing method for many different bac- 
terial species to differentiate strains and identify clonal 
lineages [6]. This technique is appropriate for global or 
long-term epidemiology and surveillance, and has the 
advantages over other typing methods, such as genetic 
fingerprinting of electronic portability and unambiguous 
characterization of isolates. MLST data are easily stored 
in databases that can be exchanged between different 
laboratories [7]. MLST has also been applied to genetic 
analysis of 160 strains of different Brucella species and 
biovars [8]. A total of 27 sequence types were defined 
among these strains. An extended MLST method was 
also developed for Brucella, with which the Brucella iso- 
lates could be differentiated with higher resolution [9]. 

To investigate etiological changes of brucellosis in 
China, in the present study, species, biovars and geno- 
types of Brucella isolates mainly from outbreaks in Inner 
Mongolia were analyzed and compared. The most com- 
mon biovars and genotypes during the three periods 
were investigated. 

Methods 

Sample collection 

All Brucella isolates were collected and preserved by the 
State Key Laboratory for Infectious Disease Prevention 
and Control, National Institute for Communicable Disease 
Control and Prevention. Brucella field strains were iso- 
lated from clinical samples from human Brucellosis out- 
breaks. The samples were collected as part of standard 



patient care between 1955 and 2008 and were fully de- 
identified. 

Bacterial strains 

All isolates were identified as Brucella species and biovar 
on the basis of classical identification procedures: C02 
requirement, H2S production, inhibition of growth by 
basic fuchsin and thionin, agglutination with monospe- 
cific antisera and phage typing. All the isolates used for 
the present study were preserved without extensive la- 
boratory passage. Brucella strains were grown to station- 
ary phase at 37°C in Tryptic Soy Broth (TSB) or Tryptic 
Soy Agar (TSA) with or without 5% CO2. Genomic 
DNA was isolated from Brucella cultures by using a 
Wizard SV Genomic DNA purification kit. 

PCR amplification and sequencing 

Nine distinct genomic loci, including seven housekeep- 
ing genes, one outer membrane protein, and one inter- 
genic fragment, were selected as described previously 
[8]. PCR reaction mixes were prepared by mixing 5 pi 
lOx PCR Buffer with MgCl 2 , 4 pi 2.5 mM dNTPs, 1 pi 
of each primer (20 pM), 0.25 pi of Pyrobest DNA poly- 
merase (5 U/pl), 33.75 pi of distilled water and 5 pi of 
diluted genomic DNA (10 ng). Cycling parameters were 
as follows: 94°C for 5 min, followed by 30 cycles of 94°C 
for 0.5 min, 63°C for 0.5 min, and 72"C for 1 min, and at 
last an elongation step of 72°C for 10 min. Products 
were separated by agarose gel electrophoresis to check 
for success of amplification and to ensure that only a 
single product of the expected size was present. PCR 
products were then purified and sequenced with ABI 
3730 with Big Dye terminator cycle sequencing kit ver- 
sion 3.1 (Applied Biosystems) at the Beijing Genomics 
Institute. 

Sequence analysis 

The raw sequence data were edited using Edit module 
and contigs were generated from forword and reverse 
sequences by using SeqMan module of the Lasergene 
package (version 5) [10]. MLST sequences (accession 
number from AM694191 through AM695630) of the 
strains described by Adrian M Whatmore were down- 
loaded from GeneBank database. Each distinct allele at 
each of the nine loci was given a numerical designation 
according to sequence of defined alleles. If the sequence 
is different from those defined previously, a new allele 
was defined. Each unique allelic pattern over all loci was 
identified as a sequence type (ST). Allelic profiles and 
sequence data were imported into the START package 
to determine mean GC content [11]. The same package 
was used to calculate the average frequencies of syn- 
onymous substitutions per potential synonymous site 
(ds) and nonsynonymous substitutions per potential 
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nonsynonymous site (cIn) by the method of Nei and 
Gojobori in order to test the degree of selection on a locus. 
A representative strain of each genotype (ST) was used for 
phylogenetic analysis. Sequences of the nine loci were then 
concatenated and used to infer phylogentic relationships 
of the representative strains. Phylogenetic analyses were 
conducted in with MEGA 5 [12]. The evolutionary rela- 
tionships were inferred using the Neighbor- Joining method 
[13]. A bootstrap consensus tree was generated to assess 
the phylogenetic support of the branching [14]. Split de- 
composition analysis of allelic profile data was performed 
using a web-based version of the SplitsTree program [15]. 

Results 

Different prominent species and biovars during three 
incidence stages 

Human brucellosis has been systematically surveyed in 
China since it was confirmed in the 1950s. According to 
the epidemiological data obtained from these surveys, 
brucellosis was highly prevalent in the 1950s and 1960s, 
and then decreased in the 1970s and 1980s period. How- 
ever, since the 1990s, the incidence increased sharply, 
forming clear three stages: high incidence (1950-1960s), 
decline (1970-1980s) and re-emergence (1990-2000s) [4]. 
Therefore, the species and biovars of the isolates were 
analyzed and compared between the three stages. A total 
of 60 strains were isolated in Inner Mongolia from 1955 
to 2008, including 24 B. abortus, 34 B. melitensis and 2 
B. suis. In total, 18, 21 and 21 strains were isolated in 
1950-1960s, 1970-1980s and 1990-2000s stages, respect- 
ively (Table 1). Species and biovars of these isolates were 
analyzed and compared according to the three incidence 
stages. As shown in Table 2, in 1950-1060 s, 3 B. abortus 
(1 biovar 1 and 3 biovar 3), 13 B. melitensis (6 biovar 2 
and 7 biovar 3) and 2 B. suis (1 biovar 1 and 1 biovar 3) 
were isolated, and in 1970-1980s, 17 B. abortus (4 biovar 
1, 11 biovar 3 and 2 biovar 6) and 4 B. melitensis (4 bio- 
var 1) were isolated, while in 1990-200 s, 4 B. abortus 
(4 biovar 3) and 17 B. melitensis (17 biovar 1) were iso- 
lated. In the high incidence stages of 1950-1960s and 
1990-2000s, the most common species was B. melitensis. 
In the low incidence stage of 1970-1980s, the most com- 
mon species was B. abortus. Furthermore, the main bio- 
vars differ between the two high incidence stages, with 
B. melitensis biovar 2 and biovar 3 in 1950-1960s and 
B. melitensis biovar 1 in 1990-2000s respectively. Dis- 
tribution difference analysis by Chi-square test showed 
that species were significantly differentially distributed 
among the three stages (j 2 = 26.425, p < 0.001). 

Different sequence type distributions among the three 
incidence stages 

To analyze the genetic polymorphism of these isolates 
by MLST, nine distinct genomic loci of the 60 strains 



were sequenced and their sequences analyzed. Each dis- 
tinct allele at each locus was given a numerical designa- 
tion and each unique allelic pattern over all nine loci 
was defined as a sequence type (ST) according to a pre- 
vious study [8]. A total of 14 distinct STs were defined. 
The most common sequence type is ST8 (24 strains, 
40.0%), followed by ST2 (10 strains, 16.7%), ST34 (9 
strains, 15.0%), ST29 (5 strains, 8.3%) and ST1 (3 strains, 
5.0%). The other 9 STs were represented by only 1 strain 
each (Table 3). ST8 is the most common sequence type 
in 1950-1960s and 1990-2000s, accounting up to 61.1% 
and 47.6% of the isolates of the stages. However, ST8 
was observed in only 14.29% of the isolates in 1970- 
1980s. Analysis with Chi-square showed that the distri- 
bution difference of ST8 among the three stages was sig- 
nificant (^ = 9.636, p = 0.006). In 1970-1980s, the most 
represented sequence type was ST2 (6 strains, 28.57%). 
The number of sequence type differed among the three 
stages, but without significance (j 2 = 3.314, p = 0.247). 
Eighteen strains from the stage of 1950-1960s belonged 
to 6 STs (3 strains per ST), and 21 strains from 1990- 
2000s belonged to 6 STs (3.5 strains per ST). However, 
the 21 strains from 1970-1980s belonged to 10 STs (2.1 
strains per ST). When analyzed in general, these STs 
were differentially distributed among the three period 
with significance (j 2 = 37.390, p = 0.007). Compared with 
sequence types of strains from other countries, 9 of the 
14 sequence types (from ST28 to ST36) were first ob- 
served by this study. This implied that sequence types of 
strain from Inner Mongolia were greatly different from 
those from other countries. 

Sequence polymorphism and alleles changes of the nine 
loci for predominant species during the three incidence 
stages 

The polymorphic alleles and sites were analyzed for the 
nine genes. The number of alleles for nine genes ranged 
from 2 (dnaK) up to 6 (aroA). The number of poly- 
morphic sites ranged from 1 {dnaK) up to 4 (omp25). 
The dN/dS ratio was calculated to determine the degree 
of selection in the sequence population. The dN/dS ratio 
ranged from 0 to 0.7678. The dN/dS ratio of dnaK and 
omp25 was 0. GIK showed the highest dN/dS ratio of 
0.761. The GC% content for the eight genes ranges from 
55.8% to 62.5%, with an average of 59.4%, being higher 
than average (57.0%) of the whole genome sequence 
(Table 4). 

To further analyze the polymorphism of the predom- 
inant species, allele changes in nine loci of the sequence 
type were compared (Table 5). In 1950-1960s, the most 
represented species was B. melitensis, and the sequence 
type is ST8 and ST34. Compared with ST8, only one site 
was changed in ST34. Allele 8 of omp25 was changed to 
allele 12. In 1990-2000s stages, the sequence types are 
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Table 1 Species, biovar, sequence type, clonal complex 
and isolation time of the 60 isolates from Inner Mongolia 

Strains ST a Time Period Stage b Species Biovar c 



Table 1 Species, biovar, sequence type, clonal complex 
and isolation time of the 60 isolates from Inner Mongolia 

(Continued) 



Bru044 


8 


1996 


1990s 


990-2000S 


B. melitensis 


1 


Bru045 


8 


1997 


1990s 


990-2000S 


B. melitensis 


' 


Bru046 


8 


1998 


1990s 


990-2000S 


B. melitensis 


1 


Bru047 


8 


1996 


1990s 


990-2000S 


B. melitensis 


1 


Bru048 


8 


1999 


1990s 


990-2000S 


B. melitensis 


1 


Bru049 


34 


2008 


2000s 


990-2000S 


B. melitensis 




Bru050 


8 


1998 


1990s 


990-2000S 


B. melitensis 




Bru051 


8 


2002 


2000s 


990-2000S 


B. melitensis 


1 


Bru052 


34 


2002 


2000s 


990-2000S 


B. melitensis 




Bru053 


35 


1997 


1990s 


990-2000S 


B. melitensis 




Bru054 


34 


2007 


2000s 


990-2000S 


B. melitensis 




Bru055 


34 


1996 


1990s 


990-2000S 


B. melitensis 




Bru056 


8 


2007 


2000s 


990-2000S 


B. melitensis 




Bru057 


8 


2006 


2000s 


990-2000S 


B. melitensis 




Bru058 


2 


1995 


1990s 


990-2000S 


B. abortus 


3 


Bru059 


29 


1996 


1990s 


990-2000S 


B. abortus 


3 


BruOeO 


2 


1997 


1990s 


990-2000S 


B. abortus 


3 



a ST represent sequence type. b Stage represent the incidence stages: high 
incidence of 1950-1 960s, low incidence stage 1970-1 980s, and re-emerging 
high incidence stage of 1990-2000s. c The Arabic number represents the 
biovar of the species. 



still most represented by ST8 and ST34, but a new se- 
quence type ST35 was observed. Compared with ST34, 
the allele of gyrB was changed from 1 to 8. However, 
only 1 strain of ST35 was observed. B. abortus is the 
predominant species at low incidence stages. In 1950- 
1960s, only 3 B. abortus strains were isolated, they 
belonged to 2 different STs, ST2 and ST5. 2 of the 9 loci 
showed allele changes between the 2 STs. Alleles of glK 
were changed from 2 to 1, and trpE from 3 to 4. In 
1970-1980s, 8 STs were observed for B. abortus. The 
most represented sequence type was ST2 (6 strains) and 
ST29 (4 strains). The two sequence types had only 1 al- 
lele differences. The allele of omp25 was changed from 1 
in ST2 to 11 in ST9. For all the 8 STs observed in 
B. abortus strains from 1970-1980s, 6 of the 9 loci chan- 
ged. A total of 19 alleles for all the loci were observed, 
indicating the high diversity of B. abortus in this inci- 
dence stages. Compared with alleles of the nine loci 
from strains of other countries, new alleles were ob- 
served in four loci, including aroA (allele 8, 9 and 10), 
cobQ (allele 8, 9), gyrB (allele 7), omp25 (allele 11, 12). 

Combined comparison of sequence type and 
species-biovar changes during the three incidence stages 

In 1950-1960s, 2 biovars of B. melitensis, biovar 2 and 3, 
were identified. The 2 biovars belonged to 2 sequence 
types, ST8 and ST34. However, there was not a strict 



Bru001 


8 


1955 


1950s 


1 950-1 960s 


8. melitensis 


2 


Bru002 


8 


1958 


1950s 


1950-1 960s 


8. melitensis 


2 


Bru003 


8 


1960 


1960s 


1950-1 960s 


8. melitensis 


2 


Bru004 


34 


1963 


1960s 


1950-1 960s 


8. melitensis 


2 


Bru005 


8 


1965 


1960s 


1950-1 960s 


8. melitensis 


2 


Bru006 


34 


1966 


1960s 


1950-1 960s 


8. melitensis 


2 


Bru007 


8 


1957 


1950s 


1950-1 960s 


8. melitensis 


3 


Bru008 


8 


1957 


1950s 


1950-1 960s 


8. melitensis 


3 


Bru009 


8 


1957 


1950s 


1950-1 960s 


8. melitensis 


3 


Bru010 


8 


1958 


1950s 


1950-1 960s 


8. melitensis 


3 


Bru01 1 


8 


1958 


1950s 


1950-1 960s 


8. melitensis 


3 


Bru012 


8 


1960 


1960s 


1950-1 960s 


8. melitensis 


3 


Bru013 


8 


1965 


1960s 


1950-1 960s 


8. melitensis 


3 


Bru014 


2 


1960 


1960s 


1950-1 960s 


B. abortus 


1 


Bru018 


5 


1956 


1950s 


1950-1 960s 


B. abortus 


3 


Bru019 


2 


1956 


1950s 


1950-1 960s 


8. abortus 


3 


Bru022 


36 


1957 


1950s 


1950-1 960s 


B. suis 


1 


Bru023 


1/ 


1963 


1960s 


1950-1 960s 


B. suis 


3 


Bru015 


28 


1972 


1970s 


1970-1 980s 


B. abortus 


1 


Bru016 


1 


1980 


1980s 


1970-1 980s 


8. abortus 


1 


Bru020 


33 


1979 


1970s 


1970-1 980s 


B. abortus 


6 


Bru021 


31 


1980 


1980s 


1970-1 980s 


B. abortus 


6 


Bru024 


2 


1988 


1980s 


1970-1 980s 


8. abortus 


3 


Bru025 


2 


1988 


1980s 


1970-1 980s 


B. abortus 


3 


Bru026 


1 


1988 


1980s 


1970-1 980s 


8. abortus 


1 


Bru027 


29 


1989 


1980s 


1970-1 980s 


B. abortus 


3 


Bru028 


29 


1985 


1980s 


1970-1 980s 


8. abortus 


3 


Bru029 


32 


1985 


1980s 


1970-1 980s 


B. abortus 


3 


Bru030 


8 


1985 


1980s 


1970-1 980s 


8. melitensis 


1 


Bru031 


2 


1985 


1980s 


1970-1 980s 


8. abortus 


3 


Bru032 


30 


1985 


1980s 


1970-1 980s 


B. abortus 


1 


Bru033 


29 


1988 


1980s 


1970-1 980s 


8. abortus 


3 


Bru034 


2 


1988 


1980s 


1970-1 980s 


B. abortus 


3 


Bru035 


2 


1988 


1980s 


1970-1 980s 


8. abortus 


3 


Bru036 


2 


1988 


1980s 


1970-1 980s 


B. abortus 


3 


Bru037 


29 


1988 


1980s 


1970-1 980s 


B. abortus 


3 


Bru038 


8 


1988 


1980s 


1970-1 980s 


8. melitensis 




Bru039 


34 


1989 


1980s 


1970-1 980s 


8. melitensis 




Bru040 


8 


1989 


1980s 


1970-1 980s 


8. melitensis 




Bru017 


1 


2006 


2000s 


1 990-2000S 


8. abortus 




Bru041 


34 


1993 


1990s 


1 990-2000S 


8. melitensis 




Bru042 


8 


1990 


1990s 


1 990-2000S 


8. melitensis 




Bru043 


34 


2001 


2000s 


1 990-2000S 


8. melitensis 
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Table 2 Species and biovar distribution during the three 
incidence stages 



Species/Biovars 


Incidence stage 3 Total 


Species Biovars 


1950-1 960s 1970-1 980s 1990-2000S 


B. melitensis 1 


0 4 17 21 


2 


6 0 0 6 


3 


7 0 0 7 


B. abortus 1 


1 4 16 


3 


2 11 3 16 


6 


0 2 0 2 


B. suis 1 


1 0 0 1 


3 


1 0 0 1 


Total 


18 21 21 60 


a The number represent the strain number of respective incidence stage. 


consistency between the STs and biovars. Both strains of 
B. melitensis biovar 2 and biovar 3 belonged to ST8 and 
ST34. For species B. abortus, 2 biovars and STs were ob- 
served. In 1970-1980s, for the most commnon species B. 
abortus, 3 biovars (1, 3 and 6) and 8 STs were observed. 
The most common sequence type and biovar were ST2 
and biovar 3. Interestingly, strains of the most repre- 
sented sequence type ST2 and ST29 belonged to biovar 
3. Therefore, in this low incidence stage, biovar 3 and se- 
quence type ST2 and ST29 are the most distributed. In 
1990-2000s, B. melitensis was the most abundant spe- 
cies, and only one biovar (biovar 1) was identified. These 
strains belonged to 3 STs, mostly represented by ST8 
and ST34, sequence type with only 1 allele different. 


Table 3 Sequence type distributions during the three 
incidence stages 


ST 


Incidence stage 3 Total 


1950-1 960s 


1970-1 980s 1990-2000S 


1 


2 (9.52%) 1 (4.76%) 3 (5.00%) 


2 2(11.11%) 


6(28.57%) 2(9.52%) 10(16.67%) 


5 1 (5.56%) 


1 (1.67%) 


8 11 (61.11%) 


3(14.29%) 10(47.62%) 24(40.00%) 


17 1 (5.56%) 


1 (1.67%) 


28 


1 (4.76%) 1 (1.67%) 


29 


4(19.05%) 1(4.76%) 5(8.33%) 


30 


1 (4.76%) 1 (1.67%) 


31 


1 (4.76%) 1 (1.67%) 


32 


1 (4.76%) 1 (1.67%) 


33 


1 (4.76%) 1 (1.67%) 


34 2 (11.11%) 


1 (4.76%) 6 (28.57%) 9 (1 5.00%) 


35 


1 (4.76%) 1 (1.67%) 


36 1 (5.56%) 


1 (1.67%) 



a The number represent strain number. The percentage in brackets represents 
the percentage of the strains in the time period. 



Compared with those in 1950-1960s, sequence type and 
biovar of B. melitensis in 1990-2000s showed an interest- 
ing change. Although still most represented by ST8, 
more ST34 strains observed in 1990-2000s. Biovar 1 in 
1990-2000s replaced biovar 2 and 3 in 1950-1960s. 

Phylogenetic relationship Brucella isolates 

To further test the phylogenetic structure of the isolates, 
sequences of the 9 loci for each of the strains were 
concatenated and analyzed in MEGA 5. All the strains 
clustered into 15 branches, representing a total of 15 se- 
quence types. As shown in Figure 1, these strains formed 
two clades. The 2 B. suis isolates form a clade. 

Discussions 

Human brucellosis in China has been documented from 
the early 1900s. The earliest two cases of brucellosis in 
China were reported to occurred in Shanghai by Boone 
in 1905 [16]. The two cases were diagnosed as having 
Malta fever by clinical symptoms and limited epidemio- 
logical data. In 1925, Y. Zongyang reported four cases 
[4]. Antibodies were detected and B. melitensis was iso- 
lated from one of the four patients. This represented the 
first confirmed brucellosis cases in China. Although sys- 
tematic surveys and effective measures were not taken 
for the prevention and control of brucellosis before 
1950, the published data suggested that both human and 
animal brucellosis were present in many provinces of 
China before 1950 [4]. Since 1950, systematic surveys of 
both human and animal brucellosis have been carried 
out. From 1952 to 1990, surveys involving about 79.7 
million animals of more than 20 species of domestic and 
wild animals showed a positive rate of 17.3% for animal 
brucellosis. Human brucellosis positive rate during this 
period was about 3.56% [17]. The general epidemio- 
logical characteristics of brucellosis in China are as fol- 
lows: in 1950-1960s, human and animal brucellosis was 
quite severe; in 1970-1980s, the incidence of human and 
animal brucellosis was relatively low, and seemed to de- 
crease during this period; in 1990-2000s, there were no 
obvious change of animal brucellosis, but the incidence 
of human brucellosis increased. Because humans are in- 
fected by Brucella through contact with infected ani- 
mals, the increase or decrease of human brucellosis rate 
usually follows the same trend of animal brucellosis [18]. 
The main reason for contrast in incidence of human and 
animal brucellosis was considered to be resulted from 
the differences in the time, place, tests and sampling 
methods of the surveys on human and animal brucel- 
losis [18]. 

Etiological analysis of human and animal brucellosis 
revealed prominent Brucella species in China [18]. Ana- 
lysis of the Brucella strains isolated from human and an- 
imals showed that the most common species was B. 
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Table 4 Sequence polymorphism of the 9 loci 



Locus 


Alleles 


Polymorphic sites 


Polymorphic sites percent 


dfg 


ds 


d^/ds 


Mean % GC 


aroA 


6 


4 


0.71% 


0.0026 


0.0058 


0.4422 


61.77 


cobQ 


5 


,\ 


0.95% 


0.0044 


0.0058 


0.7678 


59.43 


dnak 


2 


1 


0.21% 


0 


0.009 


0 


60.53 


gap 


3 


3 


0.51% 


0.0015 


0.0091 


0.1668 


58.35 


glk 


A 


4 


0.84% 


0.0029 


0.0038 


0.761 


62.53 


gyrB 


5 


3 


0.64% 


0.0017 


0.0102 


0.1681 


58.93 


omp25 


5 


4 


0.82% 


0 


0.02 


0 


59.2 


trpE 


3 


2 


0.41 % 


0.0018 


0.0056 


0.3239 


58.44 


Int-hyp 


3 


4 


0.93% 








55.81 



melitensis, followed by B. abortus and B. suis. Only a small 
proportion of the cases are caused by B. ovis and atypical 
Brucella strains [3]. B. melitensis was the most prominent 
species associated with outbreaks. But the prominent spe- 
cies differed among provinces [17]. For the regions where 
at least two species exist, the predominant strains were B. 
melitensis with either B. abortus and/or B. suis [4]. 

Inner Mongolia has experienced a high incidence in 
1950-1960s, a decrease in 1970-1980s and re-emergence 
of high incidence since 1990s [4]. Because most of the 
historical isolates were not preserved, only limited num- 
ber of strains could be included in this study. Since most 
of these strains were isolated from brucellosis outbreaks, 
it is possible to investigate general trends in terms of 
species/biovar and genotype changes. In 1950-1960s and 
1990-2000s, when the incidence was high, a higher per- 
centage of B. melitensis was observed. In 1970-1980s, 
more B. abortus were isolated. 



In 1950-1960s, 3 species and 6 biovars of Brucella were 
isolated. In 1970-1980s, 2 species {B. abortus and B. 
melitensis) and 5 biovars were isolated. In 1990-2000s, 2 
species (B. abortus and B. melitensis) but only 3 biovars 
were isolated. This indicated that in Inner Mongolia, di- 
versity in epidemic species and biovars was reduced dur- 
ing the three stages. At high incidence stages, B. melitensis 
was the most represented species. But the epidemic bio- 
vars for the two stages differed. In 1950- 1960s, 6 biovar 2 
and 7 biovar 3 strains were isolated. But in 1990-2000s, all 
17 B. melitensis strains were biovar 1. That is, the main 
biovars of B. melitensis have changed in the two time 
periods. Previous epidemiological data showed that B. 
melitensis was predominant species that was associated 
with outbreaks throughout the country [17]. Results from 
the present study further implied that different biovars of 
this species were associated with the high incidence and 
outbreaks for the two stages. 



Table 5 Allele profile and strains distribution of the 14 sequence types 


sr 


Species 










Allele distribution 13 












aroA 


cobQ 


dnak 


gap 


glk 


gyrB 


omp25 


trpE 


Int-hyp 


1 


B. abortus 






2 


2 


1 




1 


3 




2 


B. abortus 






2 


2 


2 




1 


3 




5 


B. abortus 






2 


2 


1 




1 


4 




28 


B. abortus 






2 


2 


1 




1 1 


4 




29 


B. abortus 






2 


2 


2 




1 1 


3 




30 


B. abortus 






2 


2 


3 




1 


3 




32 


B. abortus 


8 




2 


2 


2 




1 


3 




31 


B. abortus 


1 


8 


2 


2 


1 




1 


3 


2 


33 


B. abortus 


10 


9 


2 


2 


4 




1 


3 


1 


8 


B. melitensis 


2 


3 


2 


3 


3 




8 


5 


2 


34 


B. melitensis 


2 


3 


2 


3 


3 




12 


5 


2 


35 


B. melitensis 


2 


3 


2 


3 


3 


8 


12 


5 


2 


17 


B. suis 


6 


5 


1 


1 


4 


5 


2 


3 


4 


36 


B. suis 


9 


5 


1 


1 


4 


4 


2 


3 


1 


a ST represents sequence type; 


b Alleles were 


defined based 


on previous definition 
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-c 



4: 



Strain 




2 


-5 






H 


8 
1" 




int-hyp 


ST 


Species 


Biovar 


Bru014 


1.0 


1.0 


2.0 


2.0 


2.0 


1.0 


1.0 


3.0 


1.0 


2 


BA 


BA1 


Bru019 


1.0 


1.0 


2.0 


2.0 


2.0 


1.0 


1.0 


3.0 


1.0 


2 


BA 


BA3 


Bru024 


1.0 


1.0 


2.0 


2 0 


2 0 


1.0 


1 0 


3.0 


1.0 


2 


BA 


BA3 


Bru025 


1.0 


1.0 


2.0 


2.0 


2.0 


1.0 


1.0 


3.0 


1.0 


2 


BA 


BA3 


Bru031 


1.0 


10 


2 0 


2 0 


2.0 


1.0 


1.0 


3.0 


10 


2 


BA 


BA3 


Bru034 


1.0 


1.0 


2.0 


2.0 


2.0 


1.0 


1.0 


3.0 


1.0 


2 


BA 


BA3 


Bru035 


1.0 


1.0 


2.0 


2.0 


2.0 


1.0 


1.0 


3.0 


1.0 


2 


BA 


BA3 


Bru036 


10 


1.0 


2 0 


2.0 


2 0 


1.0 


1.0 


3.0 


1.0 


2 


BA 


BA3 


Bru058 


1.0 


1.0 


2.0 


2.0 


2.0 


1.0 


1.0 


3.0 


1.0 


2 


BA 


BA3 


Bru060 


1.0 


10 


2 0 


2.0 


2.0 


1.0 


1.0 


3.0 


10 


2 


BA 


BA3 


Bru029 


8 0 


1.0 


2.0 


2.0 


2.0 


1.0 


1.0 


3.0 


1.0 


32 


BA 


BA3 


Bru027 


1.0 


1.0 


2.0 


2.0 


2.0 


1.0 


11.0 


3.0 


1.0 


29 


BA 


BA3 


Bru028 


1.0 


10 


20 


2 0 


2.0 


1.0 


11 0 


3.0 


1.0 


29 


BA 


BA3 


Bru033 


1.0 


1.0 


2.0 


2.0 


2.0 


1.0 


1-.0 


3.0 


1.0 


29 


BA 


BA3 


Bru037 


1.0 


10 


2.0 


2.0 


2.0 


1.0 


11.0 


3.0 


1.0 


29 


BA 


BA3 


Bru059 


10 


1.0 


2.0 


2.0 


2.0 


1.0 


1- 0 


3.0 


1.0 


29 


BA 


BA3 


Bru016 


10 


1.0 


2.0 


2 0 


1.0 


1.0 


1.0 


3.0 


1.0 


1 


BA 


BA1 


Bru017 


1.0 


10 


2.0 


2.0 


-.0 


1.0 


1.0 


3.0 


1.0 


1 


BA 


BA1 


Bru026 


1.0 


1.0 


2.0 


2.0 


1.0 


1.0 


1.0 


3.0 


1.0 


1 


BA 


BA1 


Bru032 


10 


10 


2 0 


2 0 


3.0 


1.0 


1.0 


3.0 


1.0 


30 


BA 


BA1 


Bru021 


1.0 


s.o 


2.0 


2.0 


1.0 


1.0 


1.0 


3.0 


2.0 


31 


BA 


BAB 


Bru015 


1.0 


1.0 


2.0 


2.0 


-.0 


1.0 


11.0 


4.0 


1.0 


28 


BA 


BA1 


Bru018 


10 


1.0 


2.0 


2.0 


1.0 


1.0 


1.0 


4.0 


1.0 


5 


BA 


BA3 


A054 


3.0 


4.0 


2.0 


1.0 


9.0 


1.0 


3.0 


3.0 


1.0 


13 


B.ovis 


B.ovis 


Bru020 


10.0 


9.0 


2.0 


2.0 


4.0 


7 0 


1 0 


3 0 


10 


33 


BA 


BA6 


Bru022 


9.0 


5.0 


1.0 


1.0 


4.0 


4.0 


2.0 


3.0 


1.0 


38 


BS 


BS1 


Bru023 


6.0 


5.0 


1.0 


1.0 


4.0 


5.0 


2.0 


3.0 


4.0 


17 


BS 


BS3 


Bru004 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


12.0 


5.0 


2.0 


34 


BM 


BM2 


Bru006 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


12.0 


5.0 


2.0 


34 


BM 


BM2 


Bru039 


2 0 


3.0 


2.0 


3.0 


3 0 


1.0 


12 0 


50 


2.0 


34 


BM 


BM1 


Bru041 


2.0 


3.0 


2 0 


3 0 


3 0 


1.0 


12 0 


5.0 


2 0 


34 


BM 


BM1 


Bru043 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


12.0 


5.0 


2.0 


34 


BM 


BM1 


Bru049 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


12.0 


5 0 


2.0 


34 


BM 


BM1 


Bru052 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


12.0 


5.0 


2.0 


34 


BM 


BM1 


Bru054 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


12 0 


5.0 


2.0 


34 


BM 


BM1 


Bru055 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


12 0 


S.O 


2.0 


34 


BM 


BM1 


Bru053 


2 0 


3 0 


2.0 


3.0 


3.0 


8 0 


12 0 


5.0 


2.0 


35 


BM 


BM1 


Bru001 


2.0 


3.0 


2 0 


3.0 


3 0 


1.0 


8.0 


5.0 


2.0 


8 


BM 


BM2 


Bru002 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8.0 


5.0 


2.0 


8 


BM 


BM2 


Bru003 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8.0 


5.0 


2.0 


8 


BM 


BM2 


Bru005 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8.0 


5.0 


2.0 


8 


BM 


BM2 


BruOO? 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8.0 


5.0 


2.0 


8 


BM 


BM3 


Bru008 


2 0 


3.0 


2.0 


3.0 


3.0 


1.0 


8.0 


5.0 


2.0 


8 


BM 


BM3 


Bru009 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8.0 


5.0 


2.0 


8 


BM 


BM3 


Bru010 


2.0 


3 0 


2.0 


3.0 


3.0 


1.0 


8.0 


S.O 


2 0 


8 


BM 


BM3 


Bru011 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8.0 


5.0 


2.0 


8 


BM 


BM3 


Bru012 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8.0 


5.0 


2.0 


8 


BM 


BM3 


Bru013 


2 0 


3 0 


2 0 


3.0 


3 0 


1.0 


8.0 


5.0 


2 0 


8 


BM 


BM3 


Bru030 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8.0 


5.0 


2.0 


8 


BM 


BM1 


Bru038 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8 0 


5.0 


2.0 


8 


BM 


BM1 


Bru040 


2.0 


30 


2.0 


3.0 


3.0 


1 0 


8 0 


S.O 


20 


8 


BM 


BM1 


Bru042 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8.0 


5.0 


2.0 


8 


BM 


BM1 


Bru044 


2.0 


3.0 


2.0 


3.0 


3 0 


1.0 


so 


5.0 


2.0 


8 


BM 


BM1 


Bru045 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8.0 


5.0 


2.0 


8 


BM. 


BM1 


Bru046 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8 0 


5 0 


2.0 


8 


BM 


BM1 


Bru047 


2.0 


3.0 


2 0 


3.0 


3.0 


1.0 


8 0 


5.0 


2 0 


8 


BM 


BM1 


Bru048 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8.0 


5.0 


2.0 


8 


BM 


BM1 


Bru050 


2.0 


3.0 


2.0 


3.0 


3.0 


1 0 


8.0 


5.0 


2.0 


8 


BM 


BM1 


Bru051 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


3 0 


S.O 


2.0 


8 


BM 


BM1 


Bru056 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


8.0 


5.0 


2.0 


8 


BM 


BM1 


Bru057 


2.0 


3.0 


2.0 


3.0 


3.0 


1.0 


80 


5.0 


2.0 


8 


BM 


BM1 



Figure 1 Phylogeny based on concatenated sequence of 9 loci showing the relationships of the 60 isolates. Sequences of each sequence 
type were concatenated and phylogenetic trees were analyzed by using neighbor-joining method. 8. ovis was used as an outgroup. 
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Sequence type analysis revealed a total of 14 STs. Most 
of the STs were less common. Nine STs contained only 
one strain. The most represented ST is ST8, which con- 
tained 24 strains, accounting for 40.1% of all the tested 
strains. Comparison of ST distributions between the three 
stages showed that ST8 is the most common STs in 1950- 
1960s and 1990-2000s, two high incidence stages. ST2 is 
the most common sequence type during the low incidence 
stage 1970-1980s. All the ST2 strains belonged to species 
B. abortus. ST distribution analysis showed that 6 STs 
were observed at the two high incidence stages, but 10 
STs were observed at the low incidence stage. The reason 
for the diversity differences between sequence types of the 
strains might be the differences in species and biovars dis- 
tributions. During the high incidence stage, B. melitensis 
was the most common species, while in low incidence 
stage, B. abortus was the most represented species. Fur- 
ther analysis of the isolates showed that the number of 
intraspecies polymorphisms differed among the three spe- 
cies. B. abortus was the most genetic diverse species. A 
total of 20 polymorphic sites were observed in sequences 
of the 24 B. abortus strains. Nine STs and biovar 1, 3 and 
6 were observed for these B. abortus strains. B. melitensis 
was more homogenous, 3 STs and 3 biovars were ob- 
served. No clear relationship between biovar and ST 
within both B. melitensis and B. abortus was observed. 
This was in agreement with previous MLST results. The 
inconsistence of genotype and biovars has been observed 
in analysis with other typing approaches [2] . 

MLST is a molecular genotyping method appropriate 
for global analysis. We also compared sequence types of 
strain from Inner Mongolia and those from other coun- 
tries. Of the 14 STs, 9 (64.28%) were firstly observed in 
strains from Inner Mongolia. The ST differences were 
resulted from differences in alleles. New alleles were 
observed in 4 of the 9 loci. These data implied that 
Brucella strains from Inner Mongolia were greatly differ- 
ent from those from other countries. 

Generally, changes of genotype are not completely 
consistent with that of species or biovars of bacterial 
strains [8,19]. Combined comparison between sequence 
type and species-biovar changes during the three inci- 
dence stages showed that data from this study was also 
consistent with that. For the most represented species B. 
melitensis, in 1950-1960s, biovar 2 and 3 were mostly 
represented. These strains belonged to ST8 and ST34. 
And in 1990-2000s, biovar 1 wasmost common, the se- 
quence types were still ST8 and ST34, but more ST34 
appeared. That is, the biovars has been changed greatly, 
but the genotype remained unchanged. 

Conclusions 

In summary, in the present study, we analyzed species, 
biovars and sequence types of 60 Brucella stains isolated 



from different incidence stages in Inner Mongolia, the 
province in China with highest incidence of brucellosis. 
The species and biovars of the isolates from the high inci- 
dence stage were significantly different from that from 
low incidence stage. Genotyping by MLST of the isolates 
showed that sequence types of strains from high incidence 
stages differed from those from the low incidence stage. 
Both the biovars and genotypes in 1990-2000s were differ- 
ent from that in 1950-1960s, indicated the emergence of 
new biovars and genotypes of Brucella strains. This im- 
plied that the etiology alterations might be a reason for 
brucellosis incidence changes in China. Further epidemio- 
logical investigation will be greatly helpful for explanation 
of the etiology changes of brucellosis in China. 

Abbreviations 

MLST: Multi locus sequence typing; ST: Sequence type. 
Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

YC, YK, YW, XY, XZ, HJ performed the study and analyzed the data, ZW, QZ, 
YY, LH, BC, ZC conceived the study and analyzed the data, YC, YK, YW, ZC 
drafted and edited the manuscript. All authors read and approved the final 
manuscript. 

Acknowledgments 

This work was supported by National Basic Research Program of China 
(Grant No. 2009CB522602), National Natural Science Foundation of China 
(31000548, 31000041, 81071320, 30901071, 31272592), National Key Program 
for Infectious Diseases of China (2008ZX1 0004-01 5, 2009ZX1 0004-1 03, 
2009ZX1 0004-205, 2008ZX1 0004-008). 

Author details 

1 Department of Infectious Disease Control, Beijing Institute of Disease 
Control and Prevention, Beijing 100071, P.R. China. 2 State Key Laboratory for 
Infectious Disease Prevention and Control, National Institute for 
Communicable Disease Control and Prevention, Beijing, P.R. China. 3 School 
of Public Health, Key Laboratory of Zoonosis, Ministry of Education, Jilin 
University, Changchun, Jilin 130021, P.R. China. 4 College of Medicine, Shihezi 
University, Shihezi 832003, P.R. China. 5 Research Center of Molecular Biology, 
Inner Mongolia Medical College, Hohhot, P.R. China. 

Received: 3 August 2012 Accepted: 28 October 2013 
Published: 1 November 2013 

References 

1. Verger J, Grimont F, Grimont P, Grayon M: Brucella, a monospecific genus 
as shown by deoxyribonucleic acid hybridization. Int J Syst Evol Microbiol 
1985, 35(3)292. 

2. Moreno E, Cloeckaert A, Moriyon I: Brucella evolution and taxonomy. 
Vet Microbiol 2002, 90(1-4)209-227. 

3. Zia SH, Wang FL: Brucellosis in north China; a clinical, etiological and 
epidemiological study. Am J Trop Med Hyg 1949, 29(6)325-936. 

4. Deqiu S, Donglou X, Jiming Y: Epidemiology and control of brucellosis in 
China. Vet Microbiol 2002, 90(1 -4):1 65-1 82. 

5. Zhang WY, Guo WD, Sun SH, Jiang JF, Sun HL, Li SL, Liu W, Cao WC: Human 
brucellosis, Inner Mongolia, China. Emerg Infect Dis 2010, 16(12)2001-2003. 

6. Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE, Urwin R, Zhang Q, Zhou J, 
Zurth K, Caugant DA, et al. Multilocus sequence typing: a portable approach 
to the identification of clones within populations of pathogenic 
microorganisms. Proc Natl Acad Sci USA 1998, 95(6)3140-3145. 

7. Aanensen DM, Spratt BG: The multilocus sequence typing network: mist, 
net. Nucleic Acids Res 2005, 33(Web Server issue):W728-W733. 

8. Whatmore AM, Perrett LL, MacMillan AP: Characterisation of the genetic 
diversity of Brucella by multilocus sequencing. BMC Microbiol 2007, 7:34. 



Chen ef al. BMC Infectious Diseases 2013, 13:514 
http://www.biomedcentral.eom/1 471 -2334/1 3/5 1 4 



Page 9 of 9 



9. Chen Y, Zhen Q, Wang Y, Xu J, Sun Y, Li T, Gao L, Guo F, Wang D, Yuan X, 
ef al: Development of an extended multilocus sequence typing for 
genotyping of Brucella isolates. J Microbiol Methods 201 1, 86(2):252-254. 

10. Burland TG: DNASTAR's Lasergene sequence analysis software. Methods 
Mol Biol 2000, 132:71-91. 

1 1. Jolley K, Feil E, Chan MS, Maiden MCJ: Sequence type analysis and 
recombinational tests (START). Bioinformatics 2001, 1 7(1 2):1 230. 

1 2. Tamura K, Peterson D, Peterson N, Stecher G, Nei (VI, Kumar S: MEGA5: 
molecular evolutionary genetics analysis using maximum likelihood, 
evolutionary distance, and maximum parsimony methods. Mol Biol Evol 
2011, 28(1 0):273 1-2739. 

1 3. Saitou N, Nei (VI: The neighbor-joining method: a new method for 
reconstructing phylogenetic trees. Mol Biol Evol 1987, 4(4):406. 

14. Felsenstein J: Confidence limits on phylogenies: an approach using the 
bootstrap. Evolution 1985, 39(4):783-791. 

15. Bandelt HJ, Dress AW: Split decomposition: a new and useful approach to 
phylogenetic analysis of distance data. Mol Phylogenet Evol 1992, 

1 (3):242-252. 

1 6. Boone H: Malta fever in China. China Med Mission 1 905, 1 9:1 67-1 73. 

17. Shang D: [Progress in the study of prevention and control of Brucellosis 
in China in last 50 years]. Zhonghua liu xing bing xue za zhi =Zhonghua 
liuxingbingxue zazhi 2000, 21(1 ):55-57. 

18. Cui YL: [Current status and prospect of brucellosis research in China]. 
Zhonghua liu xing bing xue za zhi =Zhonghua liuxingbingxue zazhi 1 983, 
4(3):1 90-1 91 . 

1 9. Kattar MM, Jaafar RF, Araj GF, Le Fleche P, Matar GM, Abi Rached R, Khalife S, 
Vergnaud G: Evaluation of a multilocus variable-number tandem-repeat 
analysis scheme for typing human Brucella isolates in a region of 
brucellosis endemicity. J Clin Microbiol 2008, 46(12)3935-3940. 

• • 

doi:10.1 186/1471-2334-13-514 

Cite this article as: Chen ef al: Changes of predominant species/biovars 
and sequence types of Brucella isolates, Inner Mongolia, China. BMC 

Infectious Diseases 2013 13:514. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at /^'\ n! _, ul _-| r Q r,tr=l 

www.biomedcentral.com/submit ammBa central 



V 



